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In this paper we study research trends in condensed matter physics. Trends are analyzed by means of the 
the number of publications in the different sub-fields as function of the years. We found that many research 
topics have a similar behavior with an initial fast growth and a next slower exponential decay. We derived 
a simple model to describe this behavior and built up some predictions for future trends. 



1 Introduction 



The investigation of social phenomena has acquired a large importance in 
the last years. 1 In particular statistical physics has been applied to study 
interaction of individuals as elementary units of a social structure. In these 
models each individual interacts with a limited number of peers usually 
negligible compared with the total number of agents. In this ensemble of 
interacting agents, the statistical physics has been used to investigate opin- 
ion and cultural dynamics, crowds behavior, social spreading phenomena 
and so on. 1 

On the other hand the structure of the human connections has been sub- 
ject of an intense study by means of networks theory. 2 It has been shown 
that humans form particular kinds of networks, the so-called scale-free net- 
works, that exhibit a very short path from one site to another thanks to the 
presence of large hubs connected to many nodes. 1 Also scientists and the 
scientific community have been subject of similar studies. In particular 
network theory has been employed to show that scientific collaborations 3 
form a scale-free network. One important characteristic of these networks 
is that new ideas and trends can spread without any initial threshold due to 
the unbounded fluctuations in the connectivity. 4 

In this paper we will analyze the birth and death of trends in the scientific 
community, and in particular in condensed matter physics. In the last fifty 
years different topics have been subject of intense research in condensed 
matter physics, ranging from superconductivity, to charge density waves or 
nanostructures. Researchers, during their career, have usually worked on 
some of these subjects moving from one to another when they found the 
new one more stimulating or challenging. This collective behavior gave 
rise to the birth, growth and death of new trends. Here we want to inves- 
tigate this dynamics without dealing directly with the diffusion of ideas in 
the scientists network, but devising a simple model to describe birth and 
death of new research topics in the community. In order to do so we will 
perform a bibliometric study based on lexical analysis of papers published 
in the last fifty years. Bibliometrics has been already employed to detect 
and investigate emerging research topics 5 by mean of co-citation links or 
text-based techniques. 6 Here we will use a simple text-based approach to 
obtained the data that will then analyzed by mean of a model derived by 
population dynamics. The model resemblances the dynamics of fashion in 
social science, that has been studied by means of Langevin and Fokker- 
Planck equations. 7,8 We will show that different sub-fields, irrespective of 
their importance and time location, exhibit a similar behavior that is well 
described by our simple model. 

The questions we would like to answer are: Are researches changing topics 
more often than in the past? Do research trends last less then in the past? 
And if this is true, is it due to the larger number of researchers? We do 
not have a definite answer to all these questions but we will show that with 
a simple model we can get some insights for trends in condensed matter 
physics. The paper is organized as follows: in section[2]we discuss publi- 
cations data; in section [5] we present the model used to describe trends in 
condensed matter physics; in section|4]we apply our model to the publica- 
tions data and report the results; finally in section[5]we draw conclusions 
and discuss open questions. 
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2 Publications and data 

The number of scientific publications in physics is exponentially grow- 
ing (see Fig. \T}. This trend has been found also in other research fields 
as chemistry, biology and so on. 9 However the growth pattern of given 
sub-fields may be different from the total one, as already pointed out in 
previous studies. 9 This behavior is due to the fact that particular research 
areas are saturated, lack of interest or possible innovation. 
In order to evaluate trends in condensed matter physics we consider 
the ratio between the number of publications in a given sub-field and 
the total one. The ratio is a good indicator for the importance of a 
given research sub-field because it is independent from the total publica- 
tions number and it allows to compare research topics at different times. 
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Figure 1 : (color online) Normalization factor obtained from SAO/NASA 
Astrophysics Data System Abstract Service. 10 The total number 
of papers has been fitted with an exponential curve. Fit parame- 
ters are reported in the figure. 



We started analyzing the number of abstracts published in physics from 
1950 to 2011." We found an exponential growth with a double time of 
about 16 years (data are shown in Figure^ in agreement with previous 
studies. 12,9 Then we extracted the number of papers in the different sub- 
fields, and normalized this number with the total number of papers pub- 
lished each year. Some results are reported in Figure [3] The selection of 
papers belonging to a given sub-field has been done by means of a search 
in the abstract database using combination of several keywords to best rep- 
resent a research topic (see for instance the first column of Table[T}. 



3 The model 

In order to analyze trends in condensed matter physics we use a modi- 
fied version of the logistic equation. 13 In his original formulation Pierre- 
Francois Verhulst applied the logistic equation to the population growth, 
assuming the reproduction rate being proportional to both the existing pop- 
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ulation and the amount of available resources: 

dP ( P 

= r p 1 

dt \ K 



(1) 



where r is the growth rate and the K represents the environment carrying 
capacity. 

In our model we consider P as the fraction of scientists that decides to 
work on a given research topic and we suppose P to be proportional to the 
paper ratio. As in the logistic model we suppose that the number of sci- 
entists working on a sub-field increases due to the available resources. In 
our case the "resources" represent the possibility to make new discoveries 
in the sub-field, but also the possibility to get more grants and more cita- 
tions in an emerging research area. This gain is equivalent to the so-called 
"bandwagon" effect in fashion dynamics. 8 
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Figure 2: (color online) Simple model for fashion in physics, with parame- 
ters K(t = 0) = 10, P(t = 0) = 0.1 and parameters r = 0.1, 
s = 0.01. 

Differently from the logistic model we suppose that the carrying capac- 
ity of a given research topic decreases with time. In fact one can expect 
that after some time a large part of the possible experiments and theory will 
be already over. Moreover there is also an economical component 14 that 
contributes to the carrying capacity decrease: grants on a new and hot topic 
have clearly more successful probabilities than on an old and widely stud- 
ied subject. This fact automatically decreases the possibility to have many 
scientists working on old research topics due to the reduction of available 
positions. All these effects, that contribute to reduce the carrying capacity 
are equivalent to the so-called "snob effect" in the fashion dynamics. 
Therefore we suppose the K, the carrying capacity of our system, de- 
creases proportionally to the number of scientists working on a given re- 
search topic multiplied by a research rate coefficient s 



dK{t) 
dt 



-sP(t). 



(2) 



The research rate coefficient s represents the speed at which the properties 
K(t) are investigated by scientists P(t). In principle the decay of K (t) 
can have a more complex dependence from P(t). In literature there are 
several studies that include retardation effects or non-linear coupling, 15 
however we did not find particular reasons for using more complicated 
models than a linear one. Combining equations [T]and|2] we get: 



d 2 K(t) _ dK(t) f 1 dK(t) 1 
dt 2 ~ T dt V + ~s dt K(t) 
The solution of this equation has the form: 



P(t) 
K(t) 



' (r-s)(C 2 r+d C t ) / 



(r 
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rs 
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(3) 

(4) 
(5) 



where C\ and C2 are the integration constants that depend from P(0) 
and K(0). In order to understand the behavior of this model we consider 



the limits at short and long times of the previous solution. For t — > 00 
we found that eq[4]decays exponentially as P(t — > 00) ~ e ~ rs /( r-s )*. 
Since s is supposed to be much smaller than r the exponent reduces to 
P(t — > 00) ~ e~ st . On the other hand for short time t, P(t — > 0) goes 
as Pit — > 0) ~ c + e rt . A typical sketch of the solution is reported in the 
Figure |2] 

From the limits at short and long time we can see that the initial growth 
is dominated by the r coefficient while the exponential decay is dictated 
by the research rate coefficient s. In the language of population dynamics 
the exponential decay is related to the decrease of carrying capacity, while 
the initial growth is due to the diffusion of a research topic in the scientific 
community. In order to analyze the publication data of different sub-fields 
we simplify eq.|4]in such a way to catch short and long time behaviors with 
a simple formula: 



Pmodel(t) = Pq- 



1 



-a(i-io) + e ^(t-to) 



(6) 



This model keeps only the main exponents for the growth and the decay 
part respectively, plus a constant Prj that is the peak intensity while to is 
related to the maximum position. We added a constant c not present in the 
initial equation that will be discussed in the next section. 
Notice that this model does not apply to any research topic, in fact there 
are some active fields that continue to grow without no sign of decay, as 
for instance research on Silicon, Density Functional Theory, Nanowires 
and so on (see Figure [3). Moreover this model can be easily generalized 
to take into account more peaks as for instance in superconductivity where 
the different discoveries gave rise to a sum of multiple curves, as the one 
in Figure[3] translated in time and with different intensities. 



Topic 


to 


(P a /Psc) 


P 


a 


crystal fields 


1957.0(7) 


0.095(5) 


0.018(2) 


0.27(4) 


helium- 3/-4 


1959(2) 


0.051(3) 


0.06(1) 


0.17(5) 


alkali halides 


1960(1) 


0.186(6) 


0.070(6) 


0.12(1) 


spin waves 


1960.51(1) 


0.125(12) 


0.044(9) 


0.30(4) 


Haas-van- Alphen 


1960.9(5) 


0.023(3) 


0.060(7) 


0.5(2) 


Brullouin scattering 


1969(1) 


0.048(3) 


0.024(4) 


0.24(3) 


Compton scattering 


1971.1(3) 


0.0173(9) 


0.060(4) 


0.8(1) 


spin glass 


1977.9(3) 


0.097(3) 


0.018(2) 


0.6(1) 


charge density waves 


1978.8(7) 


0.054(3) 


0.026(4) 


0.33(4) 


Josephson junctions 


1983(1) 


0.11(1) 


0.038(8) 


0.15(1) 


tract. Quantum Hall 


1983.8(3) 


0.0162(1) 


0.018(4) 


1.6(5) 


four wave mixing 


1984.5(5) 


0.091(5) 


0.050(5) 


0.30(2) 


quasicrystals 


1985.7(2) 


0.031(2) 


0.025(4) 


1.6(5) 


Penrose Tiles 


1986.1(3) 


0.022(1) 


0.023(5) 


1.0(2) 


superconductivity 


1986.7(2) 


1.00(6) 


0.045(5) 


3(2) 


slave bosons 


1989.5(3) 


0.0133(4) 


0.062(6) 


0.66(7) 


t-J model 


1990.1(1) 


0.037(1) 


0.090(5) 


1.2(1) 


fullerens 


1991.2(1) 


0.169(4) 


0.030(2) 


2.8(4) 


Luttinger liquid 


1993.1(3) 


0.024(1) 


0.032(5) 


0.76(9) 


giant-magnetores . 


1993.0(1) 


0.049(2) 


0.071(7) 


3(1) 


Andreev reflection 


1993.0(6) 


0.013(1) 


0.014(8) 


0.54(9) 


critical exponents 


1993(2) 


0.163(7) 


0.078(1) 


0.10(1) 


Bethe ansatz 


1994(1) 


0.045(3) 


0.08(1) 


0.19(1) 


composite fermions 


1994.7(3) 


0.022(2) 


0.11(1) 


1.4(3) 


critical phenomena 


1995(2) 


0.096(4) 


0.09(2) 


0.10(1) 


spin-charge separation 


1999(1) 


0.0125(1) 


0.15(5) 


1.0(4) 


Coulomb blockade 


1999(1) 


0.069(6) 


0.08(1) 


0.26(2) 


nanotubes 


2003.1(4) 


0.57(4) 


0.08(1) 


0.41(1) 



Table 1 : (color online) Parameters obtained by fitting the normalized num- 
ber of abstracts with the model in Eq.[2] The strength of each re- 
search topic has been normalized to the one of superconductivity 
P sc = 0.0262. 



4 Data fit and results 

We suppose that the relevance of a research topic can be evaluated from 
the number of papers published on that subject. We extracted this number 
from the SAO/NASA Astrophysics Data System Abstract Service (SAS). 10 
Since the total number of papers published every year is exponentially 
growing, we normalized the number of papers of a given research area with 
the total one in order to evaluate the importance of each research topic. 
Also the normalization factor is obtained by the SAO/NASA database (see 
Fig. [T}. The data for the different topics were then fitted with the model 
described by Eq.|6] where the constant c has been always set to zero except 
for the case of "Superconductivity". In fact in this case there was already 
an active research community working on "BCS" superconductivity before 
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the discovery of high-Tc superconductors boosted again the number of pa- 
pers. In order not to miss papers and improve search results we included 
a list of synonymous for each research area as for instance "superconduc- 
tors, superconductor, superconductivity etc....". The fit results are reported 
in Table [T] where we renormalized the peak intensity with the one of "su- 
perconductivity", the larger peak in our sample, in such a way to underline 
the relevance of the each sub-field. Selected data and relative fits are also 
plotted in Fig.0 

Notice that recent research topics exhibit larger error bars because we 
could fit only the beginning of the decay tail. Moreover some active 
sub-fields as "graphene, topological insulators, quantum computing, pho- 
tovoltaic, etc..." are sill in the growing part of the curve or have just started 



to decrease and this makes impossible to perform an accurate fit with our 
model. 

In Tabled we report the different parameters entering in the model that 
describes trends in condensed matter physics. Unfortunately the results for 
the a parameters are very noisy due to the fast diffusion of new research 
topics in the scientific community. This makes more difficult to get an ac- 
curate estimate for this exponent. Nevertheless we found that a is slowing 
increasing with time(to ), and this indicates a faster spread of ideas among 
scientists. This fast spread can be explained by the fact that the structure of 
the collaboration networks is scale-free 3 and in this kind of networks there 
is a vanishing threshold to spread ideas and trends. 4 
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Figure 3: (color online) Number of papers for different research topics versus the years normalized to the total one in physics. The blue solid lines 
represent the fitting curves, fitting parameters are reported in Table[T] In green two topics that do not fit in our model. 



Then we analyze the f} parameter. In Figure[4]we report the parameter /3 
as a function of the maximum position. The maximum position is given by 
tmax = to + a ^ tog(g), however since a is usually one or two order 
of magnitude larger than /3, the maximum position can be in many cases 
approximated with to . 

The parameter /3 is related to the life-time of a research topic. We found 
that there is a variation of this parameter before and after the 1995. In 
fact we have an average life-time of about 25 years before 1995 while the 
average life-time drops to just 10 years in the last period. 



5 Conclusions and open questions 

In conclusion we present a simple model to describe trends in science 
and applied it to condensed matter physics. We found that different re- 
search topics are well represented by our model. In particular topics that 
had limited applications in research or industry. These topics show a fast 
growth followed by an exponential decay in publication ratio. On the con- 
trary other topics that became relevant for industry (as Silicon, transistors, 
lasers, etc..) and/or for basic research (as lasers, many-body perturbation 
theory, density functional theory,...) show a growth comparable or supe- 
rior to the total publications number in physics. We learn from these results 
that while new topics spread in the scientific community in few years their 
decay time is of the orders of decades. 
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Figure 4: (color online) The parameter as a function of t max that repre- 
sents the peak maximum position. 

These results could be used for example to investigate "scientific bub- 
bles" 16 or to discriminate between the importance of different sub-fields.' 7 
In fact one can expect that "bubbles" will have a different decay in time 
than real research sub-fields, and try to detect them. Then we analyzed the 
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a and (3 parameters that are related to the birth and death of a research sub- 
field. We found an increase of both (3 and a with time. Unfortunately the 
data are too noisy to model this behavior and to draw reliable conclusions. 
Therefore the question we posed in the title remains open: is research go- 
ing faster and faster? Or in another way are trends in condensed matter 
physics changing faster than in the past? We have some indications but not 
a definite answer. 

For sure there are different possibilities to improve the present results that 
we have not considered yet, as for instance to extend this research to all 
physics branches, including more sources or combining lexical analysis 
with citation metrics. We hope this paper will stimulate future research on 
diffusion, birth and depth of trends in the scientists network. 

Acknowledgments 

CA acknowledges Lorenzo Stella for useful discussions, and Elena Can- 
nuccia for carefully reading this paper. This work has not been supported 
by any funding program. 

References 

[1] Claudio Castellano, Santo Fortunato, and Vittorio Loreto. Statisti- 
cal physics of social dynamics. Rev. Mod. Phys., 81:591-646, May 
2009. 

[2] Reka Albert and Albert-Laszlo Barabasi. Statistical mechanics of 
complex networks. Rev. Mod. Phys., 74:47-97, Jan 2002. 

[3] M. E. J Newman. The structure of scientific collaboration networks. 
PNAS, 98(2):404, 2001. 

[4] Romualdo Pastor-Satorras and Alessandro Vespignani. Epidemic 
spreading in scale-free networks. Phys. Rev. Lett., 86:3200-3203, 
Apr 2001. 

[5] W. Glanzel. Bibliometric methods for detecting and analysing emerg- 
ing research topics. El profesional de la information, 21 (2): 194-201, 
2012. 

[6] M. Zitt and E. Bassecoulard. Development of a method for detection 
and trend analysis of research fronts built by lexical or cocitation 
analysis. Scientometrics, 30:333-351, 1994. 



[7] Wolfgan Weidlich. Physics and social science: the approach of syn- 
ergetics. Phys. Rep., 204:1-163, 1991. 

[8] Shoichiro Nakayama and Yasuyuki Nakamura. A fashion model with 
social interaction. Physica A: Statistical Mechanics and its Applica- 
tions, 337:625, 2004. 

[9] Jean Tague, Jamshid Beheshti, and Lorna Rees-Potter. The law of 
exponential growth: Evidence, implications and forecasts. Library 
Trends, 30:125, 1981. 

[10] ADS, CfA, and NASA. Sao/nasa astrophysics data system abstract 
service. 

[11] All data are obtained from the SAO/NASA abstract service 10 by 
mean of a Python script that automatically performs the researchs 
and parses the results. We tried also to perform research on other 
databases as ISI Web of Science, APS journals or arXiv and we 
found a similar behaviour. However data in SAO/NASA are larger 
than the other mentioned databases and this allows to reduce fluctua- 
tion among the different years. 

[12] L J Anthony, H East, and M J Slater. The growth of the literature of 
physics. Reports on Progress in Physics, 32(2):709, 1969. 

[13] Pierre-Francois Verhulst. Notice sur la loi que la population poursuit 
dans son accroissement. Correspondance mathematique et physique, 
10:113, 1838. 

[14] Stephan Paula E. The economics of science. Journal of Economic 
Literature, 34(3): 1199, 1996. 

[15] V.I. Yukalov, E.P Yukalova, and D. Sornette. Punctuated evolution 
due to delayed carrying capacity. Physica D: Nonlinear Phenomena, 
238(17): 1752- 1767, 2009. 

[16] Abraham Loeb. Rating growth of scientific knowledge and risk from 
theory bubbles, http://arxiv.org/abs/1108.5282, 2011. 

[17] Abraham Loeb. How to nurture scientific discoveries despite their 
unpredictable nature, http://arxiv.org/abs/1207.3812, 2012. 



4 



CITATION DETAILS 



